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This  piiase  of  the  Nutatron  gyroscopic  sensor  development  was  sponsored  by  the  Air  Force 
Avionics  Laboratory,  WPAFB  under  Contract  F33615-69-C-1722  and  by  Bell  Aerospace  Company. 
The  Nutatron  gyro  features  an  anisometric  rotor  which  responds  to  torques  by  exhibiting  minute 
signals  at  twice  the  rotor  frequency.  These  signals  are  used  to  automatically  compensate  for  the 
drift  producing  torques  and  as  a  result  the  Nutatron  is  a  high  performance,  low  cost  gyro  with  fast 
reaction  time.  Under  previous  contracts,  the  feasibility  of  the  concept  was  demonstrated  and  an 
engineering  model  was  designed,  fabricated,  and  tested.  This  report  covers  a  phase  of  the  instrument 
development  which  had  as  its  goal  the  investigation  of  the  effect  of  bearing  imperfections  on  the 
instrument  performance.  During  this  program  unique  analyses  of  bearing  noise  were  made  and  the 
results  of  these  analyses  show  which  bearing  parameters  are  important  to  achieving  good  Nutatron 
performance. 
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SYMBOLS  USED  IN  SECTION  II1-C 


Fa  =  Axial  Load 
Fr  =  Radial  Load 

=  Individual  Ball  Load 
a  =  Loaded  Contact  Angle 

4/  -  Ball  Position  Angle 

6a  =  Axial  displacement 

Sr  =  Radial  displacement 

kn  =  Ball  load  vs.  ball  race  deformation  coefficient 
A  =  rQ  +  r  j  -  D 

r0  =  Radius  of  curvature  of  outer  race 
rj  =  Radius  of  curvature  of  inner  race 
6  a  =  6a/A  normalized  axial  displacement 
Ir  =6  r/A  normalized  radial  displacement 
5  n  =  The  normal  approach  of  the  races 
0  =  0 /A  normalized  angular  displacement 

a0  =  free  contact  angle 

Rj  =  Radius  of  locus  of  centers  of  curvature  of  inner  race. 
R0  =  Radius  of  locus  of  centers  of  curvature  of  outer  race. 
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small  conical  radius 

large  conical  radius 

direction  of  motion  of  model  slider 

direction  perpendicular  to  motion  of  model  slider 
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axial  force 
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average  pressure  at  y,  <j> 
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SECTION  I 


INTRODUCTION 


This  document  constitutes  the  final  report  on  a  unique  bearing  study  phase  of  the  Nutatron 
development  program.  Initially  this  phase  was  to  cover  the  design  and  incorporation  of  a 
hydrodynamic  gas  bearing  in  the  Nutatron  gyro  engineering  model.  Due  to  a  change  in  the  funding 
situation,  the  work  statement  was  revised  to  cover  the  design  of  a  gas  bearing  and  a  noise  analysis  of 
the  Nutatron  ball  bearing.  The  Air  Force  Avionics  Laboratory,  sponsored  this  program  under 
Contract  No.  F33615-69-C-1722. 

The  Nutatron  is  an  angular  rate  sensor  using  an  anisometric  rotor.  The  instrument  responds  to 
gyroscopic  torques  by  exhibiting  minute  oscillations  at  twice  the  rotor  frequency  in  addition  to 
conventional  gyro  drift.  An  automatic  feedback  loop  senses  these  minute  oscillations  and  nulls  the 
drift  which  produces  them.  Hence,  continuous  drift  compensation  is  provided  and  the  requirement 
to  stabilize  torques  prior  to  a  navigation  mission  is  eliminated.  As  a  result,  the  Nutatron  is  a  high 
performance  gyro  with  an  extremely  short  readiness  time  (less  than  one  minute). 

During  previous  phases  of  the  Nutatron  development  effort,  a  Nutatron  engineering  model  was 
designed,  fabricated  and  tested.  Several  error  mechanisms  were  identified  and  corrected  so  that  the 
only  major  remaining  source  of  error  is  noise  at  twice  the  rotor  frequency  arising  from  the  ball  type 
spin  bearings.  Under  this  phase  detailed  examination  of  the  characteristics  of  the  bearing  noise  was 
undertaken.  The  steps  required  to  reduce  the  effects  of  bearing  noise  are  now  known  and  can  be 
incorporated  in  the  instrument  in  a  subsequent  effort.  Table  I  illustrates  the  effect  these  steps  will 
have.  When  these  steps  have  been  taken,  the  instrument  will  be  ready  for  an  independent  evaluation 
program  at  a  facility  such  as  Holloman  AFB. 

The  body  of  this  report  is  devoted  to  the  work  completed  under  this  phase  of  the  development 
only.  A  discussion  of  the  Nutatron  theory  and  a  description  of  the  engineering  model  can  be  found 
in  Technical  Reports  AFAL-TR-65-170,  AFAL-TR-67-54,  and  AFAL-TR-69-76. 
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SECTION  II 


SUMMARY 


The  primary  task  of  this  program  was  a  bearing  study  effort  including  the  design  of  a  low 
speed  gas  bearing  for  the  Nutatron  and  a  thorough  analysis,  supplemented  by  a  test  effort,  of  ball 
bearing  noise  specifically  at  the  Nutatron  frequency.  The  results  of  these  efforts  showed  that: 

1 .  The  spin  bearings  are  the  major  remaining  source  of  error  producing  noise.  This  noise 
affects  performance  in  four  ways; 

a)  Geometric  errors  in  the  bearing  produce  noise  at  twice  the  rotor  frequency. 
However,  this  source  has  been  shown  to  be  very  small. 

b)  Improper  mounting  of  the  bearing  is  the  major  source  of  noise  that  is  g  sensitive. 
This  can  be  corrected  with  minor  modification  to  the  rotor  structure. 

c)  Bearing  noise  occuring  around  the  spin  axis  is  cross  coupled  into  the  input  axes  by 
means  of  the  flex  leads.  The  use  of  very  compliant  flex  leads  will  eliminate  this  problem. 

d)  Beating  between  the  ball  complement  speeds  of  the  two  spin  bearings  gives  rise  to  a 
low  frequency  “random”  drift  component.  Purposely  introducing  a  difference  in  the  two  bearings 
will  raise  the  beat  frequency  so  that  it  can  be  filtered. 

2.  When  the  modifications  described  above  have  been  incorporated,  the  Nutatron  will  be 
ready  for  evaluation  by  an  independent  test  facility. 

3.  A  low  speed  gas  bearing  for  the  Nutatron  is  feasible.  Although  improved  performance  is 
anticipated,  a  thorough  study  of  gas  bearing  noise  at  the  Nutatron  frequency  should  be  undertaken 
before  such  a  bearing  is  used. 
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SECTION  III 

TECHNICAL  DISCUSSION 


The  approach  being  used  in  the  development  of  the  Nutatron  gyro  is  one  that  combines 
investigative  testing  of  an  engineering  model  supported  by  a  rigorous  theoretical  analysis  of 
potential  error  mechanisms.  When  an  error  source  is  identified  by  this  approach,  appropriate 
changes  are  made  to  either  the  instrument  or  the  electronics  to  eliminate  the  error  source. 
Calibration  and/or  trimming  to  reduce  the  error  level  are  avoided  because  these  operations  lead  to 
expensive  instruments  and  provide  a  source  of  potential  failure  if  changes  occur  due  to  aging, 
temperature,  etc. 

A.  Major  Results  of  the  Test/Analysis  Effort 

The  results  of  the  test/analysis  effort  show  that  four  sources  of  error  remain. 

1 .  Noise  at  twice  the  rotor  frequency  due  to  imperfections  in  the  bearing  races.  This  error 
source  is  not  a  limiting  factor  at  the  present  level  of  performance,  but  does  represent  an  area  of 
potential  improvement.  Analysis  of  the  bearing  errors  show  that  a  sufficient  improvement  can  be 
obtained  without  the  use  of  expensive  specially  made  bearings.  This  improvement  requires  a  change 
in  the  ball  count  and  contact  angle  of  the  spin  bearings. 

2.  Noise  at  the  Nutatron  frequency  resulting  from  improper  mounting  of  the  bearing.  This 
has  been  shown  experimentally  and  analytically  to  be  the  largest  source  of  g-sensitive  noise. 
G-sensitive  motion  at  twice  per  revolution  results  from  unequal  bearing  compliances  in  the  major 
axes  of  the  rotor  (see  Figure  1 ). 

Under  radial  acceleration,  the  bearing  deflects  along  the  acceleration  vector  by 


d  = 


Ma 

K 


where  M  =  the  rotor  mass 

a  =  acceleration 

K  =  compliance  of  bearing 

From  this  it  can  be  seen  that  if  Kx  *  Ky,  the  deflection  of  the  bearing  has  a  component  at  twice  per 

revolution  in  addition  to  a  steady-state  value.  Misalignment  of  the  bearings  results  in  unequal 
distribution  of  the  preload.  This  in  turn  produces  the  unequal  compliances  around  the  bearing  due 
to  the  inherent  nonlinearity  of  the  ball  compliance.  A  more  detailed  analysis  of  the  phenomena  is 
given  in  Section  III.C. 

An  improved  method  of  mounting  the  bearings  in  the  rotor  structure  will  result  in  a 
significant  reduction  of  this  error  source  without  increasing  the  complexity  of  the  instrument. 

3.  Cross  coupling  of  noise  about  the  spin  axis  into  the  input  axes  by  the  flex  leads.  Both 
g-insensitive  noise  and  random  noise  have  been  shown  to  arise  from  this  source.  The  cross  coupling 
effect  lias  been  traced  to  the  fact  that  the  flex  leads  do  not  lie  in  a  plane  through  the  center  of 
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Figure  J.  Unequal  Bearing  Compliance 
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suspension.  As  noise,  0^j  sin  2vt.  around  the  spin  axis  exercises  the  flex  leads  (Figure  2)  a  force  is 
developed  in  each  flex  lead  which  is 

f  =  Kp  0 n  s'n  -  vi. 


Because  the  plane  of  the  flex  leads  is  displaced  from  the  center  of  suspension  by  the  distance  C  the 
flex  lead  forces  are  converted  to  a  torques  about  an  input  axis. 

T  =  t:  Kp  0 sin  2  vt 


It  the  flex  compliances  (Kp^  Kp^  and  Kp^  are  eqUa|5  a  torque  balance  condition  exists.  Unequal 

compliances  however,  upset  this  balance  and  a  torque  at  twice  per  revolution  (2^)  is  introduced  into 
the  input  axis  producing  unwanted  noise. 

Elimination  of  this  error  source  will  be  obtained  by  modifying  the  flexure  pivot  to  permit 
the  use  of  very  compliant  flex  leads.  The  present  suspension  system  requires  stiff  flex  leads  to 
prevent  low  frequency  spin  axis  resonance. 

4.  A  beating  between  the  ball  component  frequencies  of  the  two  spin  bearings  adds  a  low 
frequency  periodicity  to  the  Nutatron  output.  Intentionally  making  the  bearings  at  each  end  of  the 
rotor  slightly  different  will  increase  the  beat  frequency  so  that  it  is  easily  filtered. 

The  subsequent  sections  of  this  report  provide  detailed  discussions  of  the  test  program  as 
well  as  the  analytical  effort. 


B.  Summary  of  the  Test  Program 

During  the  period  covered  by  the  contract,  a  continuous  investigative  test  program  was 
conducted  on  the  Nutatron  engineering  model.  The  major  effort  was  devoted  to  examining  the 
characteristics  of  the  Nutatron  noise.  This  was  accomplished  by  a  series  of  exploratory  tests  in 
conjunction  with  some  modifications  to  the  gyro  rotor.  The  parameters  tested  for  were  g-insensitive 
drift,  g-sensitive  drift,  randomness  and  the  temperature  coefficient  of  these  drifts.  The  present  level 
of  performance  of  these  parameters  is  illustrated  by  Table  I.  The  significant  results  of  the  test 
program  are  discussed  below. 

1 .  Measurement  of  g-Sensitive  Noise 

Under  normal  circumstances,  the  g-sensitive  Nutatron  noise  is  measured  by  positioning 
the  gyro  spin  axis  parallel  to  the  earth  rotation  axis.  In  this  position  approximately  0.7  g  is 
introduced  in  the  input  axis,  but  the  earth  rotation  rate  is  not  sensed.  The  gyro  is  operated  in  the 
normal  rate  constrained  mode  of  operation  with  the  case  rotating.  The  difference  between  the  noise 
level  measured  in  this  position  and  the  systematic  noise  measured  with  the  spin  axis  vertical  is  the 
g-sensitive  noise. 


5 


Figure  2.  Flex  Lead  Cross  Coupling 
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TABLE  I 

NUTATRON  PERFORMANCE  PARAMETERS 


Present 

Anticipated 

Systematic  Drift 
Repeatability 

0.016  deg/hr 

0.007  deg/hr 

Random  Drift  (1o) 

0.024  deg/hr 

<0.01  deg/hr 

g  Sensitive  Drift 

0.2  deg/hr/g 

0.05  deg/hr/g 

g  Sensitive  Drift 
Repeatability 

0.016  deg/hr/g 

0.007  deg/hr/g 

Anisoelastic  Drift 

None 

None 

Temperature  Sensitivity 

0.002  deg/hr/°C 

<0.002  deg/hr/°C 

Magnetic  Sensitivity 

None 

None 

g  Limit  Capability 

10  g 

10  g 

Reaction  Time 

<One  Minute 

<One  Minute 

Size 

3  in.  dia.  x  3  in.  1 

3  in.  dia.  x  3  in.  1 

CONFIDENTIAL 


CONFIDENTIAL 


An  alternate  method  of  measuring  the  g-sensitive  noise  has  proven  very  informative  in 
examining  the  noise  characteristics.  Again  the  gyro  is  positioned  with  the  spin  axis  parallel  to  the 
earth  rotation  vector.  In  this  test,  however,  the  instrument  case  is  not  continuously  rotated.  Instead, 
the  instrument  is  indexed  in  45°  increments  through  360°  about  the  spin  axis  and  the  amplitude 
and  phase  of  the  2v  noise  is  measured  for  each  position.  The  variations  in  the  2v  are  proportional  to 
the  g-sensitive  noise. 

Analysis  of  the  test  results  show  that: 

a)  A  sideband  signal  is  generated  that  is  proportional  to  g-sensitive  noise  but  is  not  at 
the  Nutatron  frequency. 

b)  The  sideband  signal  tracks  the  g-sensitive  noise  to  within  10%  to  20%  and  could 
eventually  be  used  as  a  rough  measure  of  g-sensitive  noise.  While  it  could  possibly  be  used  as  a 
compensation  signal,  the  direction  of  the  development  and  analysis  efforts  is  to  reduce  the 
g-sensitive  noise  to  a  level  that  requires  no  compensation. 

c)  The  g-sensitive  noise  is  not  a  simple  amplitude  modulation  of  the  2v  noise  level  by 
the  pickoffs,  suspension  system,  electronics,  etc.,  but  rather  is  noise  generated  within  the 
rotor/bearing  structure. 

d)  The  indexing  test  can  be  run  without  installing  the  rotor  in  the  case.  Therefore, 
investigation  of  g-sensitive  noise  can  be  done  with  the  rotor  in  the  open  to  permit  easy  access  to  the 
bearings,  etc. 

2.  Cross  Coupling  Test 

The  cross  coupling  of  noise  from  the  spin  axis  into  the  sensitive  axes  was  measured  by 
temporarily  adding  a  special  pickoff  to  sense  motions  around  the  spin  axis.  By  comparing  the 
output  of  this  pickoff  to  the  regular  pickoff,  the  amount  of  cross  coupling  could  be  measured.  The 
test  showed  that  much  of  the  random  noise  near  the  Nutatron  frequency  appearing  in  the  sensitive 
axis  was  cross  coupled  from  the  spin  axis.  It  was  also  determined  from  the  test  that  the  g-sensitive 
noise  originates  directly  in  the  sensitive  axes. 

3.  Effect  of  Bearing  Misalignment 

Analysis  of  the  effects  of  bearing  misalignment  showed  that  this  could  be  a  major  source 
of  g-sensitive  noise.  To  confirm  this,  a  rotor  was  disassembled  and  the  bearings  were  realigned  to 
improve  their  parallelism.  Testing  of  the  g-sensitive  noise  before  and  after  reassembly  showed  a 
significant  improvement  in  performance  (factor  of  three).  Further  improvement  could  not  be 
obtained  on  this  rotor  due  to  limitations  on  the  parallelism  achievable.  As  a  result,  the  bearing 
mounting  provisions  on  all  subsequent  rotors  will  be  modified. 

4.  Randomness  Investigation 

Close  investigation  of  the  characteristics  of  the  randomness  of  the  Nutatron  output 
revealed  a  low  frequency  periodicity.  The  frequency  of  this  signal  varies  greatly  from  rotor  to  rotor. 
The  source  of  this  signal  is  the  beating  of  the  ball  complement  frequencies  of  the  two  spin  bearings. 
The  more  closely  matched  the  bearings  are,  the  lower  the  beat  frequency.  Eventually  bearings 
dissimilar  enough  to  give  a  high  frequency  beat  will  be  used  so  this  signal  can  be  filtered. 
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5.  Temperature  Test 

(C)  The  internal  temperature  of  the  gyro  was  varied  by  20°C  and  the  various  drift  parameters 
were  measured.  The  results  showed  a  drift  coefficient  of  0.002°/hr/°C  for  all  drift  parameters.  It 
was  also  shown  that  the  bearing  alignment  affects  the  temperature  coefficient  of  g-sensitive  drift.  It 
is  anticipated  that  improved  alignment  will  reduce  this  temperature  coefficient  to  less  than 
0.001  °/hr/°C. 

6.  Independence  from  Operating  Parameters 

(U)  A  series  of  tests  were  run  to  determine  the  dependence  of  the  performance  parameters  on 
various  operating  conditions  such  as  electronic  loop  gains,  etc.  The  tests  showed  that  the  same 
performance  is  obtained  independent  of: 

a)  Pickoff  excitation  voltage 

b)  Pickoff  gaps 

c)  Constrainment  loop  response 

d)  Position  of  the  rotor  in  the  case 

e)  Level  of  mass  unbalance 

f)  Dynamic  unbalance  of  the  rotor 

g)  Housing  rotation  rate 

h)  Resonance  conditions 

i)  Axis  alignment  loop  response 


C.  Analysis  of  Ball  Bearing  Noise 

(U)  The  analysis  of  ball  bearing  noise  at  twice  the  rotor  frequency  was  focused  on  two  types  of 
geometric  errors: 

1 .  An  elliptical  runout  of  the  outer  (spinning)  bearing  race, 

2.  A  tilt  of  the  outer  race  due  to  mounting  imperfections.  The  analysis  showed  that  both  of 
these  errors  produce  g-sensitive  noise  at  twice  the  rotor  frequency.  The  noise  arising  from 
a  mounting  tilt,  however,  was  shown  to  be  much  higher  than  the  noise  due  to  race 
imperfections.  The  symbols  used  in  the  analyses  are  listed  following  the  Table  of 
Contents.  The  details  of  the  analysis  follow. 

1.  Effect  of  Elliptical  Races 

(U)  This  analysis  is  simplified  by  the  assumption  that  all  balls  remain  in  a  state  of 
compression  for  all  possible  conditions  of  loading.  This  requirement  is  satisfied  in  the  instrument  by 
setting  an  initial  bearing  preload  of  approximately  3  pounds.  Intuitively,  if  one  or  more  balls 
become  unloaded  during  a  revolution  of  the  rotor,  additional  undesirable  noise  can  be  expected  to 
appear. 

(U)  For  calculation  of  displacements,  the  mathematical  approach  used  was  to  assume  a 
displacement,  calculate  the  individual  ball  loads  and  sum  these  up  to  obtain  the  resultant  loads.  A 
digital  computer  was  used  to  perform  the  iterative  calculations. 
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a)  Derivation  of  Equations 

Starting  with  the  equations  for  the  sum  of  radial  and  axial  loads  (Reference  1): 


Fr  =  ]L  Qi/,  cos  a  cos  i p 

\p  =  — 7T 


(4) 


\j/  =  +  n 

Fa  =  H  sin  a 

0  =  ~7T  w 


(5) 


The  individual  ball  loads  are  found  in  terms  of  the  heaviest  loaded  ball. 


®\}/  ~  Qmax  £  1  +  2e  (  *  ‘  cos  'Z'  )J 


3/J 


•i  ('♦*?=■) 


(6) 


(7) 


In  order  to  evaluate  Qmax,  the  normal  ball-race  deformation,  5max,  is  required 


^max  =  ^a  sin  a  +  cos  a 
Qmax  =  Kn  6max3/2 


(8) 

(9) 


Kn  may  be  found  using  equations  contained  in  Reference  1  or  in  Reference  2.  Typically  for 
bearings  of  the  size  used  in  the  Nutatron,  Kn  varies  from  0.3xl07  to  O.SxlO7  lb/in.3  n . 

b.  Calculation  of  Displacements 

Using  the  foregoing  equations,  a  computer  program  was  written  to  calculate  the 
radial  motion  of  the  bearing  races  relative  to  each  other  for  constant  radial  load.  The  motions  were 
calculated  for  a  Ig  radial  load  and  for  zero  radial  load.  The  difference  between  these  yields  the  “g” 
sensitive  noise  attributable  to  race  ellipticity. 

Figure  3  is  a  flow  chart  illustrating  the  computer  solution  for  the  noise. 

c.  Discussion  of  Results 

Results  of  the  computer  analysis  for  four  types  of  bearings  are  tabulated  in  Table  II. 
An  elliptical  error  in  the  bearing  race  produces  a  small  g*sensitive  2v  motion.  The  2v  motion  results 
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Figure  3  Flow  Chart  for  Elliptical  Race  Solution 


1 1 


TABLE  II 

EFFECT  OF  ELLIPTICAL  GEOMETRY  ERRORS 
ON  RADIAL  DISPLACEMENT 


Item 

New 

Hampshire 

Standard 

Miniature 

Precision 

Bearing 

Standard 

New 

Hampshire 

Modified 

Miniature 

Precision 

Bearing 

Modified 

No.  of  Balls 

8 

10 

11 

13 

Ball  Size 

0.0781  in. 

0.0625  in. 

0.0781  in. 

0.0625  in. 

Free  Contact  Angle 

15.0° 

15.0° 

12.5° 

12.5° 

Typical  2  i>Noise 
(u  sec) 

705 

560 

500 

460 
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because  the  elliptical  error  produces  a  modulation  of  the  spring  rate  at  twice  per  revolution  due  to 
nonlinearity  of  the  ball  versus  race  deformation  equation.  However,  because  the  bearing  is 
preloaded,  the  percent  modulation  of  the  ball  load  is  relativ  ly  low  and  the  amplitude  of  the 
g-sensitive  motion  is  small  compared  to  the  elliptical  race  error. 

The  expected  geometry  error  to  be  encountered  in  actual  bearings  is  difficult  to 
predict.  Usually  the  “out-of*round”  errors  are  specified  and  measured  as  total  errors  without  regard 
to  the  number  of  cycles  around  the  circumference.  This  makes  it  difficult  to  separate  the  elliptical 
error  from  the  other  geometrical  errors.  Typically,  an  ABEC  7  bearing  will  be  specified  as  being 
100  microinches  out  of  round  on  the  inner  race  and  200  microinches  out  of  round  on  the  outer 
race.  The  bearing  manufacturers  contracted  (New  Hampshire  Bearings  and  Miniature  Precision 
Bearings  Inc.)  believe  that  the  out-of-round  errors  are  less  than  30  microinches  consisting  largely  of 
third  and  fifth  harmonics.  They  declined  to  comment  on  probable  amount  of  second  harmonic  as  it 
is  concealed  by  the  odd  harmonics. 

2.  Effect  of  Tilted  Races 

As  was  assumed  in  the  analysis  above,  all  balls  remain  in  a  state  of  compressive  loading. 
The  following  analysis  is  based  upon  the  same  conditions  as  for  the  elliptical  race  errors  so  that 
direct  comparison  can  be  made. 

a.  Derivation  of  Equations 

For  calculation  of  the  effect  of  a  tilted  race,  it  is  necessary  to  use  an  analysis  which 
takes  into  account  a  moment  loading  in  addition  to  the  axial  and  radial  loads.  The  basic  equations 
are  derived  in  Reference  1 . 

The  normal  approach  of  the  races.  6n.  is  modified  in  this  analysis  to  include  a  tilted 
outer  race  and  angles  of  tilt  not  in  the  plane  of  the  radial  load.  Reference  1  gives  the  following 
equation  for  the  inner  race  tilted  in  the  plane  of  the  radial  load: 


6n  =  A  j  [  (  sin  a  °  +  5a  +  Rj  0  co$4>  )2  +  (  cos  a  °  +  6r  cos  \p  )2  )  1  2  -1 }  (10) 

For  a  tilted  outer  race,  RQ  is  substituted  for  Rj  everywhere  it  appears.  (RQ  and  Rj  are  the  locus  of 
centers  of  curvature  of  the  bearing  races.)  If  the  outer  race  is  tilted  in  a  plane  other  than  the  plane 
of  the  radial  displacement,  the  equation  is  modified  by  introducing  a  phase  angle  in  the  6  cos  ^ 
tern.  The  normal  approach  equations  are  summarized  for  the  various  options: 

Inner  Race  Tilted: 

In  Plane  of  Radial  Load: 

6n  =  A  |  [  (  sin  a  0  +  6a  +  Rj  fl  cos  i//  )J  +  (  cos  a  0  +  Sr  cos  \!/)2  ]  1  2  -1 } 

In  Perpendicular  Plane: 

6n  =  A  { [  (  sin  a  0  +  6a  +  Rj  0  sin  \p  ) 2  +  (cos  a  0  +  5r  cos  \p)2  ]  - 1 }  (11) 
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Outer  Race  Tilted: 

In  Plane  of  Radial  Load 

5n  =  A  { ( ( sin  a 0  +  5a  +  RQ  0  cos  t//  ) 2  +  (  cos  a 0  +  6f  cos  ^  ) 2  ]  1 11  - 1 }  (12) 

In  Perpendicular  Plane: 

6n  =  A  |  ((  sin  a  0  +  6a  +  RQ  0  sin  \p  ) 2  +  (  cos  a  0  +  6r  cos  ^  ) 2  ]  1/2  -if 
The  appropriate  equation  can  be  used  to  calculate  the  normal  ball  load. 

Q  =  Mn  *  S 

The  individual  loaded  ball  contact  angles  may  be  found  from: 


sin  a 0  +  5a  +  Rj  0  cos  ^ 

sin  a  =  - — - — - — -  (14) 

(  (  sin  a  0  +  6a  +  Rj  0  cos  ) 2  +  (  cos  a  0  +  5r  cos  \p)*  ]  1/2 

cos  a  0  +  6r  cos  \p 

cos  a  =  - - - - . .  — .  . ■—  —  (15) 

[  ( sin  a  0  +  6y  +  Rj  0  cos  i//  ) 2  +  (  cos  a  0  +  1>T  cos  \J/ ) 2  ]  1/2 

The  appropriate  substitutions  for  Rj  and  0  cos  ^  must  be  made  as  before. 

To  find  the  resultant  radial,  axial,  and  moment  loads,  the  following  summations  must 

be  performed. 

-  2  7T 

Fa  =  sin  a  (16) 

^  =  0 
=  2ir 

Fr  =  £ Q0  cos  *  cosa  (17) 

^  =  0 
1^  =  2  TT 

M,  *  Vidj^Q^  cos  i// sin  a  (18) 

t  =  0 
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\i/  =  2  jt 


M2  =  */2  dm^Q^  sin  ^  sin  a  (19) 

$  =  0 

In  most  cases  either  M,  or  M2  will  he  equal  to  zero. 

b.  Calculation  of  Displacements 

The  method  of  calculation  was  to  find  6a  by  iteration  assuming  a  preload  Fa  and  an 
angular  tilt  0  and  setting  5f  =  0.  Then  6r  was  found  by  iteration  for  a  radial  one  “g”  load  and  5a 
constant.  A  new  Fa  and  the  moments  were  then  calculated.  This  computational  technique  is 
illustrated  by  the  flow  chart  (Figure  4). 

c.  Discussion  of  Results 

Examination  of  the  results  (Figures  5  and  6)  shows  a  definite  preponderance  of  2v 
noise  when  the  non-g-sensitive  noise  is  subtracted  out.  It  is  of  interest  to  note  that  the 
non-g-sensitive  noise  is  predominantly  at  the  spin  frequency  while  the  g-sensitive  noise  is 
predominantly  at  2i>.  As  the  noise  amplitude  is  proportional  to  the  square  of  the  angle  of  outer  race 
tilt,  it  can  be  concluded  that  improvement  of  the  alignment  of  the  races  relative  to  each  other  is  of 
primary  consideration. 


D.  Model  Equation  for  Bearing  Noise 

The  test  program  showed  that  essentially  there  are  four  types  of  bearing  noise  at  the  2v 
frequency. 

1 )  g-insensitive  coning  noise  at  2v.  This  noise  can  be  represented  by  a  vector  that  rotates  in 
the  same  direction  as  the  rotor. 

2)  g-insensitive  anticoning  noise  at  2v.  In  this  case,  the  noise  vector  rotates  counter  to  the 
rotor. 

3)  g-sensitive  2v  noise  that  cones  at  the  case  rotation  frequency  SI. 

4)  g-sensitive  2v  noise  that  anticones  at  the  case  rotation  frequency. 

For  the  purpose  of  investigating  these  noise  components,  a  bearing  race  model  equation  was 
derived.  This  model  demonstrates  the  types  of  errors  in  the  bearing  that  are  required  to  produce  the 
various  types  of  noise  mentioned  above.  This  analysis  also  verified  to  a  large  extent  the  conclusions 
drawn  from  the  other  bearing  analyses,  viz.  that  bearing  mounting  imperfections  have  a  greater 
impact  on  Nutatron  performance  than  geometric  errors  within  the  bearing. 

The  model  used  considers  the  bearing  races  to  be  distorted  on  their  radii,  and  a  general  Fourier 
series  representation  is  assumed.  The  spring  constant  between  races  is  assumed  nonlinear.  The  three 
dominant  bearing  race  errors  are  shown  to  be: 

1 )  An  elliptical  outer  race  distortion  coupled  with  a  first  harmonic  inner  race  distortion. 

2)  An  elliptical  outer  race  distortion  coupled  with  a  third  harmonic  inner  race  distortion. 

3)  An  elliptical  outer  race  distortion  coupled  with  a  g-sensitive  displacement  of  outer  rotor 
with  respect  to  inner  rotor. 
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TABLE  III 

PARAMETERS  USED  IN  THE  ANALYSIS 


Bearing  Data: 

Mean  Diameter  =  Dm  =  0.375  in. 

Ball  Diameter  =  D  =  0.0781  in. 

No.  of  Balls  =  Z  =  8 

Free  Contact  Angle  =  a°  =  15.5° 

Race  Curvature  Ratio  =  fG  =  fj  =  0.555 

Ball/Race  Spring  Rate  Coefficient  =  Km  =  0.304  x  107  lb/in.3 

Loads: 

Preload  (Axial)  =  Fa  =  3  lb 

Radial  Load  at  1  g  =  0.165  lb 

Angular  Tilts  =  0.0001  to  0.0005  radian 
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g  Sensitive  Motion, 
Microarcseconds 


Rad  Tilt  of  Outer  Race 


Figure  5  g  -  Sensitive  2v  Motions  versus  Angle  of  Rotation 
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Amplitude  of  g  Sensitive  Effects 
Microarcseconds 


The  same  race  geometry  errors  give  rise  to  both  straight  2v  noise  and  g-modulated  7v  noise. 

Only  one  bearing  race  geometry  is  considered;  the  other  spin  bearing  is  assumed  ideal.  The 
bearing  outer  race  translational  motions  give  rise  to  angular  motions  and  the  total  angular  motions 
are  a  superposition  of  the  motions  of  both  bearing  outer  races  assuming  the  opposing  bearing  in 
each  case  to  be  ideal. 

1 .  Discussion  of  the  Results  of  the  Analysis  (Appendix  A). 

The  bearing  induced  angular  2v  noise  (from  one  bearing)  can  be  expressed  as: 


'nx 


r0j 
R  K, 


^  sin  (2  Pt  +  )  +  6  sin  (2  vt  +  fit  +  *  -2  02 ) 

H2WK  1 
1  +  cos  fit! 


(20) 


ffk2  r0 


ny 


2  o,  r  it 

R  K  I r'3  cos  *  -  +  Pi )  “  5  cos  (2  id  +  f2t  +  ^  ~2  (32 ) 


+  r;  cos  ( 2  vt  + 

li 


*■>]  [ 


2WKj 
1  +  ■  sin  ftt 
K,2 


1 


(21) 


where: 


®nx'  %y 

R 

K,  = 
K2  = 
k2 
W 


5 

fit 


Pi 

Pi 

Hi 

a,  = 

<*3 


angular  motions  between  races  in  the  x  and  y  axes,  respectively. 

distance  between  bearings  (cm) 

linear  translation  spring  constant  of  bearing  (dyne/cm) 

Second  order  translational  spring  constant  of  bearing  (dyne/cm2 ) 

Second  order  radial  pressure  constant  of  bearing  (dyne/cm2 ) 

Rotor  weight  (dyne) 

First  and  third  harmonic  distortion  of  inner  race  (cm) 

Second  harmonic  distortion  of  outer  race  (cm) 

W 

Translational  displacement  of  outer  race  due  to  rotor  weight  (6  =  p-  )  (cm) 

Ki 

Angle  of  g-vector  with  respect  to  an  inner  race  coordinate  system.  (The  inner 
race  is  considered  stationary  and  the  g*vector  rotates  opposite  to  direction  of 
case  rotation), 
a,  -  20j 
3a3  -  202 

Angle  of  second  harmonic  outer  race  distortion 
Angle  of  first  harmonic  inner  race  distortion 
Angle  of  third  harmonic  inner  race  distortion 
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The  bearing  induced  noise  consists  of  a  coning  and  an  anticoning  component;  the  terms 
multiplied  by  r^  and  5  in  equations  (20)  and  (21 )  are  coning  at  2i>,  whereas  the  terms  multiplied  by 
r^  are  anticoning  at  2v. 

The  pickoff  response  to  these  noise  inputs  arc: 


V 

=  A 

ari, 

0, 

y_ 

f  sin  (2  v\  +  0,  -tj)*  I"  -  sin  (2  vt  +  02 ) " 

i,  +  B  ri 

'  L~COS  (2  Ul  +  T?)  J  3  _  “  COS  (2  l>t  +  0j  ). 


fib, 


/  (" sin  (2  in  +  0,  -  tj) 

alr»i  cos  0 

1  r  sin  (2  Pt  +  02  -  TJ)-] 

+  ri3  sin  0 

\  L-COS  (2  Vt  +  0,  -  TJ 

y  L-C°s(2  v\  +02  -rjy 

B|  I"  sill  (2  vl  +  0  +  03f| 
yi  B:  [--cos  (2  vt  +  0  +  03  )J 


(22) 


+  sin  0 


where: 

2jrki  r., 

A  =  - - 

R  K, 

2mk2  W  r 

B|  =  - - 

R  K,2 

2trk,  Wrc  K, 

B,  =  - — 

R  K,3 

_ _ 

J(K~2v2  A  )2  +4v2  Rd2 

v2  (2C  +  A')(l  +  N) 
b  =  - 

K  -  2  y2  (2C  +  A) (I  +P) 

A'  =  A  +  B  -  C 
A  =  A  +  B-C  +  21 


v  =  operating  speed 

K  =  pivot  spring  constant 
A,B,C  =  rotor  inertias 

I  =  non-rotating  inertia 

Rj  =  damping  constant 

2  V  Ri 

tj  =  tan'1  - 

K-2  v2  A 

P,  N  =  rotor  geometry  constants 

0  =  nt-f 

03  =  —  -2  p2 
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The  terms  multiplied  by  A  are  non*g-sensitive  2v  responses.  The  term  multiplied  by  rj 
cones  at  2v,  whereas  the  term  multiplied  by  anticones  at  2v.  The  terms  multiplied  by  y/2  B2  are 
g-sensitive  2v  responses.  Of  these,  the  terms  multiplied  by  the  constant,  a,  cone  at  2v  and  the  terms 
multiplied  by  b  anticone  at  2v. 

As  the  case  is  rotated  through  360°,  fit  moves  from  0  to  360°  or  0  rotates  from  -45°  to 
315°.  The  resultant  ellipses  in  the  x  and  y  axes  when  plotted  in  polar  coordinates  either  cone  or 
anticone  depending  on  the  relative  magnitudes  of  r^,  r^  and  BJy/2  B2,  and  on  the  relationship 
between  the  phase  angles,  ,  y2 ,  <^3  and  rj.  For  example,  let  us  assume  that  B|  =  0  =  b,  then  the 
pick  off  responses  reduce  to: 


[sin  (2  vt  +  -  r?) 

-COS  (2  Vt+ifit  -  tj) 


+  y/2  B2  a  |  cos 


[sin  (2  vt  +  v?,  -  i?)"l 
-cos  (2  vt  +  -  r?)J 


T  sin  (2  vt  +  <^2  -  tj) 
+  r:  sin0  I 

3  L-cos(2  vt  +^j  -1?) 


(23) 


The  first  term  is  a  resonant  non-g-sensitive  2v  noise.  The  second  term  represents  a  resonant  coning 
or  anticoning  ellipse  when  these  responses  are  demodulated  at  2v.  If  1 px  and  <p2  are  in  the  same  half 
plane  « 90°  apart),  the  resultant  ellipse  cones  with  respect  to  the  case  angle,  f2t  =  0  +  7r/4. 
Conversely,  if  <^,  and  *p2  are  *n  opposite  half  planes  (>  90°  apart)  the  resultant  ellipse  anticones.  If 
ri,  =  ri3  an(*  =  ^2>  the  result  is  a  coning  circle,  and  if  t£,  =1 p2  +  tt,  the  result  is  an  anticoning 
circle. 


In  the  above  example  we  have  assumed  that  the  term  multiplied  by  Bi  is  negligible.  This 
term  is  a  direct  g-sensitive  noise  which  cones  both  at  2v  and  fl.  The  ratio  Bi/^2  B2  =  K.Jy/2  K2 
=  rQ  has  the  dimensions  of  length.  When  this  term  is  included  with  the  other  terms  multiplied  by  the 
resonant  factor,  a,  it  is  difficult  to  predict  the  necessary  relative  magntiudes  of  r^ ,  rj3  and  rQ,  and 
the  relative  phase  angles  $ , ,  <fi2  and  >fi3  for  coning  or  anticoning  at  SI. 

The  terms  multiplied  by  b  are  nonresonant,  anticoning  (at  2v)  responses.  For  the 
non-g-sensitive  part  of  the  response  (terms  multiplied  by  A  in  Equation  (22) ),  the  part  of  the 
response  multiplied  by  b  gives  non-g-scnsitivc  vectors  (after  demodulation  at  2v)  which  are  not  at 
right  angles  in  the  x  and  y  axes.  The  tangent  of  the  angle  between  these  vectors  can  be  expressed  as: 


tan  (  6y  -  0X  )  = 


A,2  ~  A; 2 

2A,  Aj  sin  ($2  -1/),  +  77) 


(24) 
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where 


0X,0^  =  polar  angles  of  x  and  y  axis  vectors 

A,  =  An  a 

M 

A:  =  A  r.  b 

*3 

6 y  -  0X  =  90°  only  if  either  A2  =  0  or  =  ^2  +  *?•  If  A2  =£  0,  there  exists  an  anticoning  component 
at  2v.  In  this  model,  the  anticoning  component  is  proportional  to  rj  and  is  not  resonance 
amplified,  whereas,  the  coning  component  is  proportional  to  r^  and  is  resonance  amplified.  The 
non  perpendicularity  of  the  two  non-g-sensitive  noise  vectors  depends  on  the  relation  between  A, 
and  A2  or  in  part  on  the  relation  between  the  nonresonant  response  gain,  b  and  the  resonant  gain,  a. 

The  inner  race  distortions  can  possibly  be  caused  by  a  bump  or  a  series  of  bumps.  Assume 
the  inner  race  distortion  is  caused  by  one  bump  which  can  be  represented  by  the  peak  of  a  sine 
wave  as  illustrated  in  Figure  7.  If  the  circle  is  opened  as  shown  in  Figure  8,  a  Fourier  series  can  be 
written  for  the  harmonic  content  of  the  bump: 


r 


7r  ( 1  -  cos4-)  [ 


''sin  (n+1  )■£■ 


sin  ( n—  1 )  4 

— iFT - 


.  n  <j)  n  0 
-  2  sin  -*r-  cos  tt- 


) 


cos  n  0  + 


(25) 


where 

0  =  segment  angle  of  bump 

A  =  height  of  bump 

0  =  angle  of  fundamental  wave  period 

We  are  interested  in  the  amplitudes  of  the  fundamental  and  third  harmonic  of  the 
fundamental. 


1st: 


4  “  sin  4  cos  4 


TT  (  I  -  COS  ) 


3rd: 


•  bill  -r 


COS 


3 


TT  COS 


<P 

T 
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Assume  0  to  be  small  enough  such  that  cos  0  = 
amplitudes  are: 


1  and  sin  0  =  0  -  ,  then  the  Fourier 


1st 


4 -($-£)  ('-40 


“8 


03  A  8 
12ir  02 


1  M 

3  IT 


2  0-§£  ^ 


3rd: 


b  2  /30  27  3\  /  0*\  A  ^  03A  8  _  2  A0 
:  -I  --45  3- 

T 


for  small  0,  the  fundamental  and  third  harmonic  Fourier  amplitudes  are  equal  and  in  phase. 


In  a  practical  bearing,  one  would  not  expect  a  single  bump  on  the  inner  race  but  a 
continuous  distortion  around  the  perimeter  of  the  race.  Based  on  the  analysis  of  the  single  bump, 
however,  one  would  expect  the  amplitudes  of  the  fundamental  and  third  harmonic  distortion  to  be 
within  the  same  order  of  magnitude. 


The  validity  of  this  analysis  depends  on  the  following  beating  errors: 


( 1 )  fundamental  geometry  errors  on  inner  race  (bump) 

(2)  third  harmonic  geometry  errors  on  inner  race 

(3)  second  harmonic  geometry  errors  of  outer  race 

(4)  g-sensitive  relative  translation  of  races 

(5 )  a  second-order  bearing  translational  spring  constant. 

There  is  no  question  concerning  the  existence  of  the  first  four.  From  theoretical  bearing  studies, 
there  is  a  question  regarding  the  existence  of  a  second-order  bearing  translational  spring  constant  of 
the  order  necessary  to  give  the  2v  noise  levels  observed  during  the  test  program.  The  second  order 
spring  constant  necessary  to  produce  performance  errors  of  the  magnitude  of  those  observed  in  test 
is  100  times  what  is  felt  to  be  a  reasonable  value.  Therefore,  it  can  be  concluded  that  the  noise  due 
to  geometric  errors  in  the  bearing  races  is  at  least  one  order  of  magnitude,  and  very  nearly  two 
orders,  below  the  present  performance  level.  Furthermore,  the  use  of  special  made  ultrahigh 
precision  bearings  is  unnecessary  because  the  precision  of  conventional  instrument  bearings  is 
sufficient  for  a  high  performance  Nutatron. 


E.  Gas  Bearing  Study 

The  long  life  and  low  random  noise  level  inherent  in  a  gas  lubricated  bearing  made  the 
application  of  a  gas  bearing  to  the  Nutatron  look  attractive.  For  this  reason,  the  initial  work 
statement  for  this  contract  included  the  design,  fabrication  and  test  of  a  low  speed  hydrodynamic 
gas  bearing  for  the  Nutatron  rotor.  Due  to  a  change  in  the  funding  status  of  the  contract,  only  the 
preliminary  design  effort  was  completed. 

1.  Methods  of  Analysis 

The  origi  ■!  method  of  analysis  of  the  gas  bearings  was  adapted  from  the  thrust  bearing 
analysis  contained  ...  References  3  and  4.  Later  developments  in  the  field  of  gas  bearings  permitted 
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(U)the  use  of  the  design  charts  contained  in  Reference  5.  As  can  be  seen  from  Table  IV.5,  the 
simplified  approach  is  adequate  for  most  practical  design  calculations. 

(U)  A  disadvantage  of  the  design  charts  is  that  the  geometries  are  optimized  for  maximum 
stiffness  and  for  isoelastic  rotor  support.  This  limits  their  use  for  optimizing  2v  noise  levels  at  the 
expense  of  isoelasticity. 

2.  Gas  Bearing  Design 

(U)  The  preferred  configuration  for  a  gas  lubricated  spin  bearing  is  a  co-apex  opposed  conical 
configuration.  This  configuration  provides  a  design  in  which  the  air  gap  is  insensitive  to  temperature 
differences  and  in  general  minor  geometrical  variations  will  affect  both  radial  and  axial  loads  in  the 
same  manner.  The  spiral  groove  type  of  bearing  was  selected  because  the  average  support  pressures 
are  higher  than  for  other  types.  These  higher  support  pressures  yield  higher  load  capacities,  lower 
noise  levels  and  smaller  radial  and  axial  displacements. 

(C)  The  design  parameters  contained  in  Table  IV  were  obtained  by  using  the  following  initial 
figures: 


Ambient  Pressure 
Equivalent  Gap 
Outside  Radius 
Inside  Radius 
Gas  Viscosity 
Rotor  speed 


1.02  x  106  dyne/cm2 
1.27  x  10-4  cm 
1.22  cm 
0.648  cm 

1.4  x  I0"4  dyne  sec/cm2 
105  rad/sec 


These  figures  were  based  on  using  the  largest  bearing  geometry  that  would  fit  in  the  present 
Nutatron  housing,  the  minimum  practical  running  clearance  (50/tiin)  and  the  slowest  practical 
running  speed  ( 1 000  rpm ). 
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TABLE  IV 

NUTATRON  GAS  BEARING  PARAMETERS 


Parameter 

Calculated 

BAC  No.  60003-142 

From 

MTI  68TR29 

1.  Spiral  Groove  Angle 

75° 

72.3° 

2.  Groove  Depth  Ratio 

2.0 

2.06 

3.  Axial  Spring 

Constant 

2.6  x  1010  dynes/cm 

1.43  x  1 0* 0  dynes/cm 

4.  Groove  Length  Ratio 

0.7 

0.7 

5.  Attitude  Angle 

N.A. 

o 

CM 

6.  Anisoelastic  Apex 

Angle 

35.3° 

32.5° 

7.  Groove  Depth 

2.54  x  10"4  cm 

2.62  x  10-4  cm 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  analyses  and  investigative  tests  conducted  under  this  program  have  produced  the  following 
significant  conclusions  and  recommendations: 

A.  Conclusions 

1.  Expensive  special  made  bearings  are  unnecessary  because  noise  due  to  imperfections  in 
the  ball  bearing  races  is  not  limiting  the  instrument  performance. 

2.  G-sensitive  Nutatron  noise  results  from  improper  bearing  alignment. 

3.  Proper  bearing  alignment  can  be  achieved  without  increasing  the  instrument  complexity. 

4.  Cross  coupling  of  noise  from  the  spin  axis  into  the  input  axes  can  be  eliminated  by  using 
more  compliant  flex  leads. 

5.  By  using  two  spin  bearings  with  different  ball  complement  speeds,  the  beat  frequency 
component  of  the  randomness  can  be  eliminated. 

6.  After  modifying  the  engineering  model  as  indicated  above,  it  will  be  ready  for  test  at 
Holloman  AFB  (or  any  other  suitable  facility). 

7.  A  low  speed  gas  bearing  is  feasible. 

8.  A  thorough  noise  study  is  necessary  before  a  gas  bearing  is  incorporated. 


B  Recommendations 

* 

1.  Modify  the  Nutatron  rotor  structure  to  insure  proper  bearing  alignment. 

2.  Incorporate  flex  leads  that  are  very  compliant. 

3.  Use  spin  bearings  with  slightly  different  contact  angles  to  eliminate  beat  frequency 
problem. 

4.  Evaluate  the  modified  engineering  model  Nutatron  at  an  independent  test  facility  such  as 
HoJIoman  AFB. 

5.  Incorporation  of  a  low  speed  gas  bearing  should  be  postponed  pending  further 
investigation. 
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APPENDIX  A 


Analysis  of  a  Bearing  Model  Equation 


Referring  to  Figure  9,  let  the  inner  race  be  stationary  with  an  x-y  coordinate  system  tied  to  it. 
The  outer  race  rotates  CCW  at  v  and  let  the  g-vector  rotate  clockwise  at  ft.  (Corresponds  to  a  CCW 
case  rotation.) 

Let  r,  =  r,o  +  r,(  cos  (0  +  a, )  +  r^  cos  2  (0  +  a2 )  +  r^  cos  3  (0  +  a3)  +  . . . 

ro  “  ri„  +  d  +  ro,  cos  (0O  +  >3, )  +  r0j  cos  2(0O  +  02 )  +  r0j  cos  3  (0Q  +  p3 )  +  . . . 

+  5  cos  (0  -  \l/) 


where 


OQ 

6 

W 

0 

v 

n 


nominal  inner  race  radius 
ball  diameter 

j1*1  harmonic  distortion  of  inner  race  radius 
j1*1  harmonic  distortion  of  outer  race  radius 
j1*1  distortion  phase  angle  of  inner  race 
j1*1  distortion  phase  angle  of  outer  race 
outer  rotor  displacement  from  rotor  weight 
rotor  weight 
angle  of  g-vector 

angle  variable  in  stationary  coordinates 
angular  velocity  of  outer  race 
angular  velocity  of  g-vector 


The  outer  race  is  moving  CCW,  then  the  instantaneous  gap  between  races  as  a  function  of  8  is 


Ar  =  r0  -  Tj  and  =  0  -  v t 

Assume  a  linear  continuous  pressure  distribution  radially  outward  on  the  outer  race  according 
to  the  equation: 

Pr  =  Po+Pi(tj  rQ  +  d)  + pj  (rrr0 +  d)2  +p3  (rrr0  +  d)3  +  ...  (26) 


Then  the  incremental  radial  outward  force  is: 
dfr  =  Pr  roo  60 


where 


Let 


oo 

Po  = 

Pj  = 

S  = 


r,0  +  d  is  the  nominal  outer  race  radius 

radial  preload  pressure 

j1^1  order  of  radial  pressure 

ro0Pj  be  the  jdl  order  of  radial  spring  constant. 


Figure  9  Sketch  of  Geometry  of  Inner  and  Outer  Bearing  Race 


Then  expanding 


dfr  =  jk0  +k,  |rj(  cos  (0  +  a, )  -  rQ(  cos  (0  +  /3,  -  vt) 

+  cos  2  (0  +  a2 )  -  rQi  cos  2  (0  +  p2  -  vt) 

+  cos  3  (0  +  a3)  -  rQj  cos  3  (0  +  03  -  vt)  +  .  . . 

+  6  cos  (0-^)1  (27) 

+  k2  (  r(i  cos  (0  +  a, )  -  r0|  cos  (0  +  /3,  -  vt) 

+  cos  2  ( 0  +a2 )  -  r0j  cos  2  (0  +  ji2  -  vt) 

+  rj  cos  3  (0  +  a3 )  -  ra  cos  3  (0  +  03  -  vt)  +  . . . 

+  5  cos  ( 0  —  i^)| 2 
+  k3  I  ...l3  +  ...  }d0 

We  are  interested  only  in  force  terms  which  give  rise  to  translational  motions  of  the  outer  race 
with  respect  to  the  inner.  The  net  forces  along  the  x  and  y  axes  are  found  by  intergrating  on  0  from 
0  to  2rr  as  follows: 


fx  = 

f2ir  dfr 

l  M 

cos  0  d  0  +  W  cos  ^ 

fy  = 

f2*  dtr 
1  d0 

sin  0  d  0  +  W  sin  ^ 

where 


C/  =  -  (fit  +  tr/2)* 

fx  =  x-axis  force 

fy  =  y-axis  force 


Since  dff/d0  is  multiplied  by  either  cos  0  or  sin  0  in  equation  (28),  the  integrals  are  non-zero  only 
when  the  terms  in  dff/d0  are  also  multiplied  by  either  cos  0  or  sin  0.  Further,  only  terms  containing 
2vt  are  to  be  considered.  The  first  few  of  the  terms  meeting  these  requirements  are: 


—I  =  <-2  ki  [r:  cos  ( 0  +  a , )  r,„  cos  2  (0  +  02  -  vt) 

d0  '  *  '  J 

+  r^  cos  3  (0  +  a2 )  r0j  cos  2  (0  +  02  -  vt) 

-  6  cos  (0  +  fi  t  +4  )  rGj  cos  2  (0  +  02  -  vt)  | 

+  3  k3  (r^  cos  (0  +  a,  )  rQi  2  cos2  (0  +  (3 ,  -  vt) 

+  Tj  cos  3  (0  +  a3 )  r0(  2  cos2  (0  +  (3,  -  vt) 

+  6  cos  (0  -  ^)  r0(  2  cos2  (0  +  /3,  -  vt)  ]  j 

•At  t  =  0,  the  g-vcctor  is  assumed  downward.  The  g-vector  rotates  clockwise. 


(29) 


31 


If  the  geometry  radius  errors  arc  assumed  small,  only  the  terms  multiplied  by  k2  need  be 
considered.  Expanding  these  and  ignoring  terms  involving  3  9  we  have 


dfr'"  =  ”k2  lri(  rQi  cos  (0  +  2  -a,  -  2  t>t) 

d0 

+  r(j  r0j  cos  (0  +  3  a,  -  2  02  +  2  vi) 

-  6  rGj  cos  (0  +  2  0,  -n  t  -  2  »>t)  ) 

=  _k2 1  (ri)  r0j  cos  (2  »»t  +  sp, )  +  ri}  r0j  cos  (2  in  +  <f2 ) 

-6  r0j  cos  ( 2  vi  +  fit  +  %  -  2  02 )  1  cos  0 

+  [ri(  r0j  sin  (2  vt  +  )  -  r^  rGj  sin  (2  i>t  +  ^>2 ) 

-  6  r0j  sin  (2t't  +  fit  +  *~2/32)  ]  sin 

where 


=  a.  “  2/3j 

=  3  a3  -  20, 


since 


/2  IT  ,  flit 

cos3  0  d  0  =  j  sin3  0  d  0  =  jr 

O  ■'o 


f2  * 

and  /  cos  0  sin  0  d  0  =  0 
Jo 


(30) 


*x  ~  ~  k2  ro2lri,  cos  <-  vl  +  Vi  I  +  ri3  cos  <2  vl  +  ifj ) 

-  5  cos  (2  i>t  +  fit  -  2/?2)  |  +  W  cos  (fit  +  *  )  (31) 

fy  =  -rrk2r0j  [r^  sin  (2  vt  +  *>, )  -  rj}  sin  (2  v\  +  <p2 ) 

-  6  sin  (2  »>t  +  fit  -  2  02 )  !  -Wsin(fit+5) 


The  outer  race  motions  caused  by  the  forces  in  equation  (31 )  are  determined  by  the  effective 
spring  constant  of  the  bearing.  (Inertial  and  damping  forces  are  considered  negligible  compared  to 
the  spring  forces.)  Assume  the  motion  of  the  outer  race  to  be  defined  by  the  equations 


fx  =  K,  6X  +  K2  6x3  +  . .  . 

ty  =  K,  5y  +  K2  5y3  +  .  .  . 


(32) 
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If  higher  order  than  the  2nd  are  neglected,  we  can  solve  equation  (32)  for  6X  and  6 


K,3 


5 


y 


$ 


t/Jvj 

K,3 


(33) 


The  first  terms  on  the  right  hand  side  of  equation  (33)  give  rise  to  non-g-sensitive  2v  motions  and 
the  2nd  terms  give  rise  to  g-sensitive  2v  motions  by  virtue  of  the  cross  product  terms  (W-terms) 
when  fx  and  fy  are  squared.  Including  only  the  dominant  terms  at  2v  and  2v  modulated  at  SI  we 
have 


- 


^r1(  cos  (2  vl  +  )  +  rj3  cos  (2  vl  +  ^2) 

■y/  2  W  K2  \ 

-  5  cos  (2  nt  +  lit  +  -  2  + - ‘ —  sin  fttj 

'2  r°2  l"rj  sin  (2  v\  +  )  -  r;  sin(2H+^2) 

C,  L  1  3 

/  2WK  \ 

-  8  sin  (  2  vl  +  Six  +  *  -  2  M  +  — cos  fitj 


-jr  kj  r 

K. 


(34) 


The  relationships  between  angular  motions  and  the  linear  motions  are 

6y 

^nx  R 

6x 

n  =  ~— 

Mny  R 


(35) 


where  R  =  radius  arm  to  opposing  bearing.  The  bearing  noise  angles  are 


jt  k2  r0,  f 

I  ria  sin  (2  vt +^2)-rij  sin  (2  vt  +  if, ) 

/  2  W  K  \ 

+  6  sin  (2  *»t  +  flt  +  *  -  2  /3j)J  ^1  + - -  cos  ftt J 


(36a) 


it  k2  r0,  r 

'  R  jr*  k  cos  (2  vl  +  <f2)  +  cos  (2i^t  +  if, ) 
-  8  cos  (2  vl  +  fit  +  *  -  2  02)j  ^1  + — — j-  sin 


(3bb) 


These  are  equations  ( 1 )  and  (2)  in  Section  II. 

The  pickoff  responses  to  the  inputs  can  be  found  easiest  by  first  transforming  the  noise  angles 
into  rotating  rotor  coordinates.  Considering  the  non-g-sensitive  noise  first,  the  bearing  noise  in 
rotating  coordinates  is 


°n'x  fC0SI't  sini,tl  Kx 


0ny  ^1-^1  cosidj  0ny  ^ 


A  cos  3  vl  sin  3  vl  "  cos  (  *  -  if  2) 

—  —  j*.  4m  / 

2  13  -sin  3  vl  cos  3t>t  sin  (*-f2) 


fcosi’t  -sinid"!  cos  (<f ,  +  *) 


sin  vl  cos  i>tj  sin  (f ,  + 


where  A  = 


2ir  k2  r0j 


We  inay  write  the  pickoff  response  as 


=  A  I  a 


sin  ( 2  i>t  + 

""’I 

L-cos  ( 2  vl  + 


*1 

*2  “t?)J 


[sin  (2  »>t  +  if2 ) 
cos  (2  vt  +  if2 ) 


where 


_ v 2  A  ' _ _ 

f(  K  -  2  v2  A  )2  +  4  v2  RjJ 


TT  -  the  resonant  gain 


u  _  V2  (2C  +  A  )  (  I  +  N)  .u 
b  -  - ,  the  nonresonant  gam 

K-2r>:  (2C  +  A)  (I  +P) 
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These  are  the  non-g-sensitive  parts  of  Equation  (22)  Section  II. D. 


Considering  next  the  g-sensitive  terms  in  rotating  coordinates,  we  have 


"0nxl  J  cos,,t  sin,;tl  pnxl 

0ny  2^±n  'rsinvt  cos  l,t-‘  IvJ 


2t»±n 


V2B: 


+  r; 


sin  (3  m  +  <fii)  sin  (fit  +  +  sin  (pt  +  i^j )  sin  (fit  + 

cos  (3  vt  +  ^j)sin(flt  +  -cos  (i>t  +^2)  sin  (fit 

-sin  (3  i>t  +  ifi, )  sin  (fit  +~)  -sin (vt  +ifii )  sin  (fit +  ^ 
-cos  (3  i/t  +  \fix )  sin  (fit  +  +  cos  (t>t  +  )  sin  (fit  +  ^ 


Let 

0  = 

fit  - 

Then 

sin  ( 

fit  + 

V 

B, 

+  ~ 


sin  ^t  +  fit  +•7-  2 
-cos  (i>t  +  fit  +  *  -  2  p2j 


0 1 


nx 


O’r 


ny 


+  rj  sin  0 


p 

/ 

/  4  7T\ 

\ 

/ 

CC.f^j  +- 

cos(*,  +  “  ) 

\ 

I  r.  sin  0 

V  -  / 

+  r:  COS  0 

l  3 

1 1 

/  .  IT  \ 

I 

\ 

sin  (*2  +  T J 

s*n  (^1  +2) 

/ 

\ 

“j 

• 

B, 


+  “ 


cos  i>t  -  sin  i>t 

cos(0  +  J  -  2  p2) 

sin  i>t  cos  vt 

sin  p  +  J  “  2  P2j 

(39) 
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The  pickoff  response  is 


V 

/ 

sin  (2  vl  +  ifii  -  17) 

=  a 

n/J  Bj  (  rj(  cos  <(> 

0V 

2v±Sl  ! 

-cos  (2  vt  +  \fii  -  17) 

l  vj 

\ 

where 


(40) 


These  are  the  g-sensitive  portions  of  Equation  (22)  Section  II. D. 

The  angle  17  is  included  in  the  trigonometric  terms  which  are  multiplied  by  the  constant,  a. 
This  constant  resonates  depending  on  the  spring  compensation.  If  the  system  were  undamped, 
a  would  assume  a  positive  or  negative  value  depending  on  the  magnitude  of  2v2  A  relative  to  K.  Near 
resonance,  a  approaches  -90°  in  phase.  To  account  for  the  phase  shift  in  a,  we  can  use  the  results  of 
a  previous  analysis  on  Nutatron  damping.  These  results  show  that  a  can  he  considered  a  positive 
constant  plus  a  lagging  phase  angle  included  in  all  associated  trigonometric  terms.  For  small  pivot 
spring  constants,  the  phase  of  a  approaches  -180°;  therefore,  r?— ►  1 80° .  At  resonance,  17  =  90°,  and 
for  large  K,  rj-K)°.  The  constant,  h.  is  also  modified  somewhat  by  the  damping  constant,  R.  hut  this 
term  is  far  from  any  resonance;  hence  h  <  0  can  he  calculated  with  sufficient  accuracy  independent 
of  R. 


In  practice,  the  pickoff  signals  are  demodulated  at  2v  using  both  a  sine  and  cosine  phase.  The 
phase  of  the  demodulator  references  is  arbitrary  and  unknown,  hut  the  x  axis  is  assumed  to  he  the 
East-West  axis  and  the  reference  arbitrary  phase  is  assumed  to  be  0°  with  respect  to  this  coordinate 
system  when  time  starts.  The  introduction  of  a  reference  phase  would  not  alter  the  relationship 
between  the  non-g-sensitive  and  g-sensitive  portions  of  the  noise. 

The  race  distortion  angles,  a,.  a3,  are  therefore  defined  with  respect  to  a  coordinate 
system  whose  x  axis  is  the  axis  of  zero  reference  phase. 

Demodulating  the  non-g-sensitive  portion  of  the  pickoff  signals  first. 
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(41) 


<>x  = 

A(  sin  (ifii  -tj)  -  A2  sin  (^2 

at  cos  phase 

= 

At  cos(^|  -  tj)  -  A2  cos  <fi2 

at  sin  phase 

“At  COS(i0,  —  Tj)  —  A2  COS  (^2 

at  cos  phase 

= 

At  sin  (if,  -tj)  +A2  sin  1^2 

at  sin  phase 

where 

A,  =  A  rj(  a 
A,  =  Ar:  b 

*  *3 


(42) 


The  angle  between  the  two  demodulated  vectors  0X'  and  0y'  is  found  by  first  finding  the  angle 
of  each  vector  in  the  reference  coordinate  system  followed  by  finding  the  tangent  of  the  angle 
between  them.  Let 


tan'1 

A| 

cos  x  -  Aj  cos  ifi3 

Ai 

sin  x  -  A2  sin  ifi2 

tan-1 

Ai 

sin  x  +  A2  sin  tfi2 

“At 

cos  x  -  A2  COS  ifi2 

where  x  =  -  tj 


Then 


tan  (02  -  0, ) 


tan  0 2  -  tan  0, 

1  +  tan  0j  tan  0, 


At  sin  x  +  \2  sin  _  At  cos  x  -  \2  cos  *p2 
-A,  cos  x -Aj  cos  At  sin  x  -  A2  sin  i p2 

1  +(  A|  sin  x  +  Aa  sin  \p2  \  /At  cos  x  -  \2  cos  i p2 

\-At  cos  x  -A:  cos  /  \  At  sin  x  -  A2  sin  <p2 

After  clearing  of  fractions  and  combining  terms 


A, 2 


-A, 


2 


tan  (02 


2A,  A 2  sin  {\fi2  ~  x) 


Demodulating  the  g-sensitive  portion  of  the  pickoff  signals,  we  have  in  the  x  axis: 
=  j  j^a  (rit  sin(^,-rj)-r0sin(^,-r?)) -bri3  sin^jj  cos0 

+  ^a  (rl3  sin  (^jj-17)  -  rQ  cos  (1^3-r?))  -  b  r^  sin^,J  sin  0  j  at  cos  phase 

=  \fl  Bj  j  £J  (ri,  cos  <^1  ~»?)  “  r0  cos  “  b  rj3  cos  ifit  cos  0 

+  £a  ^r,3  cos  +  rQ  sin  (^3-17))  -  b  rjj  cos  sin  0  j  at  sin  phase 


where  r_  =  ~z — 


V2K: 


and  has  dimensions  of  length. 


Similarly,  in  the  y  axis  of  the  Nutatron  we  can  write 


0”  =  y/l  Bj  |  [  a  ( -  r,t  cos  (^,  -r?)  +  rQ  cos  (^3-r?))  -b  r^  cos  02J  cos  0  (45) 

+  |a  (-  rj3  cos  (^j-17)  -  rQ  sin  (^3-1?))  -  b  r^  cos  sin  0  j  at  cos  phase 

VTBj  |  [^a  (r^  sin  _r(j  sin  (^3-17))  +  b  sin^jj  cos0  (46) 

+  £a  ( r^  sin  (^j-r?)  ~  r0  cos  (^3-rj)  )  +  b  r^  sin  J  sin  0  j  at  sin  phase 

Equations  (43).  (44),  (45),  and  (46)  are  parametric  equations  of  ellipses  in  the  x  and  y  axes. 
The  angle  variable  is  0.  Another  set  of  parametric  equations  in  terms  of  the  semi-major  and 
semi-minor  axes  of  an  ellipse  are 

("xl  _  rcosT  “S»n  ->1  f"E  cos  (Sit  +  5)1  ^ 

[_y  J  |_sin  7  cos  7J  [_±F  sin  (fit  +  5)J 


where 

x.  y  =  x  and  y  axis  coordinates  in  reference  system 

E  =  semi-major  axis  of  ellipse 

F  =  semi-minor  axis  of  ellipse 

Tit  =  0  +  rr/4 

7  =  angle  of  ellipse  semi-major  axis  with  respect  to  the  reference  coordinates 

8  =  phase  angle  of  Sit  at  points  where  the  semi-major  axis  intersects  the  ellipse. 

The  7  matrix  rotates  the  ellipse  coordinates  into  the  reference  coordinates. 

By  manipulating  these  equations,  we  can  write  x  and  y  in  terms  of  the  parameter  angle,  0.  The 
results  are 
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(48) 


M  j"  G,  cos  (7  ±  e, )  Gj  cos(7  iejfir  cos^l 
|_yj  I  LG,  sin(7±C|)  Gj  sin  (7  ±  e,  )J  0  J 


where 


Tej+fj  Ej-Fj  .  _1 

-  — 3 —  ~  ~r~  s,n  28  I 

-  sjn:5j 

1  [r  tan  1 4  +  5)] 


1  >2 


1  /  2 


tan 


tan' 1 


[if  ™«<5  *  «>] 


The  upper  signs  in  equations  (47)  and  (48)  represent  a  coning  ellipse,  and  the  lower  signs  an 
anti-coning  ellipse. 

If  the  angle  6s  0,  then  e,  =  e2  =  tan'1  F/E  and  G|  =  G2  =[(EJ+F2V2]|/J .  We  can  now  make 
identifications  between  equation  (48)  and  equations  (43)  and  (44)  in  the  x  axis  and  between 
equation  (48)  and  equations  (45)  and  (46)  in  the  y  axis  for  the  g-sensitive  portion  only.  Then 


(a) 

y/1  B2  (a  r(i  sin  (*,  -rj)  -  a  rQ  sin  (*3-r?)  -  b  rj}  sin  *2  ) 

- 

Gx 

cos (7X  ± ex) 

(b) 

y/l  B2  (arj3  sin^j-T^-arQCOs^-TjJ-br^  sin  ) 

= 

-°x 

COS (7X  T  6X) 

(c) 

VTb2  (a  r(j  cos  (^,-rj)  -  a  rQ  cos  (<^3  — r?)  -  b  rj}  cos*2  ) 

l  = 

cx 

sin  (7X  ±  ex) 

(d) 

\ZlB2(ari3  cos(s(J2-T?)  +  ar0sin  (^j-rjJ-br^  cos^,  ) 

= 

-°x 

sin  (7X  *fx) 

(e) 

v/J B2  (-ar,(  cos(^,-7j)  +  ar0cos(^3-r?)-  b  r^  cos^j) 

1  = 

Gy 

COS  (7y  1  fy) 

(f) 

B2  ( -a  r,  cosl^j-rjJ-a^sinf^s-^J-br^  cos  ) 

= 

-°y 

cos (7y  T  ey) 

<g) 

B2  (a  rjj  sin  (*,  -r?)  -  a  r0  sin  (*,-??)  +  b  r^  sin  *2  ) 

= 

°y 

sin  (7y  ±  ey) 

(h) 

\F~  B2  (artj  sin(^2-r?)-ar0cos(^3-rj)+  brj|  sin  <fit  ) 

= 

-Gy 

sin  (7y  T  ey) 

(49) 


where  the  subscripts  x  and  y  apply  to  the  x  or  y  axis,  respectively.  The  right  hand  sides  of 
equation  (49)  are  known  from  plots  of  experiments.  Let  these  knowns  be  designated  C, ,  C2,  C3, 
C4.  C5 .  C6,  (,,  C g,  then  the  procedure  for  evaluating  the  unknowns  in  these  8  equations  is  as 
follows 

Consider  the  difference  of  equations  (a)  and  (g)  and  the  sum  of  equations  (c)  and  (e) 
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(i)  -2  y/l  Bj  br,3  sin  ^  =  C,  -C, 

(j)  -2  V3*Bj  b  rj3  cos^j  =  Cj+Cs 

Taking  the  ratio  j4-| 

c,-c7 

tan  =  Ci+£}  (50) 

Similarly,  consider  equations  (b)  -  (h)  and  (d)  +  ( 0 

(k)  —2  \/3  B2  b  rtj  sin  sPi  =  Cj-Cg 

(l)  -2  y/2  B2  b  rj(  cos  =  C4+C6 

Then 

C2  -C, 

,J"*'  =  I^Tc,  <5I) 


Further  consider  equations  (a) + (g)  and  (c)  -  (e) 


(m)  2\/3  Bj  a(rj(  sin  (*,  -r?)  -  r(>  sin  (*>j-t?) ) 

=  C,  +C7 

A 

(n)  2v/2B2a(ri|  cos  (*,  -rj)  -  r()  cos  <*3-r?)) 

=  C,-C, 

and  equations  (b)  +  (h)  and  (d)  -  < f ) 

(o)  2>/2Bj  a(r,3  sin  (*j -r?)  -  rG  cos  l^3-r?)) 

=  Cj+C8 

Ip)  2\/2  B2  a(r,3  cos(^j-rj)  +  r0sin(^j  rj)  )  =  C4  -  C6 


Fliminating  r0  terms  from  these  four  equations  we  have  after  adding  equations  On)  and  (p)  and 
subtracting  equations  (n)  and  (o). 


(q)  rj( 

sin  x  +  rj  cos  y  =  S, 

(r)  r,( 

cos  x  -  rj3  sin  y  =  S} 

where 

- 

X 

=  ^i  -»?  s, 

(C,  +  c7  +  C4  -  C6  ) 

:V2Bj  a 

y 

=  ~V 

Sj 

C,  -C,  -Cj  -cB 

2\^Bj  a 
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After  some  algebraic  manipulation  we  can  solve  equations  (q)  and  (r)  simultaneously  for  x  and  y  in 
terms  of  S, ,  S2 ,  r;  and  r;  .  The  results  are 

1 1  *3 


(s)  cos  (y  + 


S5  - 


TS  r 


+  r; 


»3 


2 


(t)  sin  (x  +  i^) 
where 


SJ  +  ri, 2 
2Sri, 


S  =ys, 5  +SjJ 


Let 


=  tan'1  Sj/S, 


t  =  S/r;  z  =  y  + 

>3 

Tj 

r  =  — -  w  =  x  +  y 

r>3 

Then  from  equations  (s)  and  (t) 

(u)  2t  cos  z  =  t5  ~  r5  +  1 
< v )  2t  sin  w  =  t5  +  r5  +  1 

Solving  (u)  and  (v)  simultaneously  cv  eliminating  t  we  have 
(w)  cos5  z  =  I  -  r5  cos5  w 


Using  the  trigonometric  identity :  cos5  z  =  Vi  +  Vi  cos  2z 
and  letting  w  =  z+4  we  have, 

r5  -  I 

sin  (2z  -  5)  =  -  (52) 

V  I  +  2  r5  cos  X  +  r4 

where 

t  .  i  I  +  r5  cos  X 
o  =  tan  1  - 

r5  sin  X 

A  knowledge  of  r  and  X  allows  us  to  calculate  z  and  w,  and  knowledge  of  and  *fi2  from  equations 
(50)  and  (51)  allows  a  calculation  of  rj  and  t. 
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Finally,  rearranging  equations  (o)  and  (p),  we  have 


-<c,  +C,) 

(o)  rQ  cos  (^3  “  rj)  *  ^  ^  —  +  rij  sin  y 

Ca-C6 

(P)  rQ  sin  (.^j  -  rj >  =  ^  "  ri}Cosy 


Taking  the  ratio 


1£> 

(O) 


Cg-C, 


tan  (sP3  — 1?> 


2y/2ti 


-  cos  y 


a  ari. 


(53) 


Cj+Cg 


2y/l  B 


a  a  ri. 


■  +  sin  y 


(xt  t  = 


S  [<C,  +  C,  +  C4  -  C6  )2  +  (C3  -  C,  -  C:  -  Cg  )2] 


1  ■: 


■|3  -V-  °2  a  r 


2>/2  B- 


A  knowledge  of  t  allows  computation  of  <f3  -tj.  Then  a  knowledge  of  rj  provides  us  with  =  3W4 
from  which  dj  is  calculated.  Finally  a,  =  «£,  +  2^2  and  ofj  =  ^2+2dj/3  completing  the 
evaluation  of  the  angular  constants. 

There  remains  the  calculation  of  distortion  magnitudes  and  spring  constants  of  the  hearing. 
From  equation  (o) 


.  ,  'o 
to)  — 

'3 


where 


-(C:  +Cg) 

2yfl  B,  an 


+  sin  y 


cos  <*3  -  17 ) 


.  .  o 

,y,  - 


B, 


K, 


r: 


>3 


n/2Kj  r: 


and  is  known 


troin  equation  ( x) 

(x)  2y/2  Bj  a  r(j 

where  e  and  t  are  known  at  this  point.  Substituting  (y )  in  (x) 


c  *^\/3  k2  rp;  K ;  W  a  rj3 
^  '  R  K,3 


(/) 


4  ff  k3  r0?  K}  W  a  rtj  ^ 
RK.'l,  ‘ 


Equations  (41 )  and  (42)  allow  us  to  calculate  the  constants  A,  and  Aa  for  the  non-g-sensitive 
noise.  Then 


(aa)  -s—  =  A  r: 
b 


Substituting  (aa)  into  (z) 
2  Wa  A, 

,bb)  bKi  (r0,'rij)rij 


2»k*  roa  ri3 
RK, 


Solving  (bb)  for  r,  we  have 

2Wa  Ajt 
b  c  K,  (r0/rij) 


(54) 


We  have  a  single  equation  in  two  unknowns,  viz,  r^  and  K, .  The  value  for  K,  can  be  supplied 
independently  from  a  bearing  analysis.  The  value  for  the  Nutatron  bearing  is 

K,  =  4.63  x  I010  dyne/cm. 


Now  r.  can  be  calculated. 


=  K,/v/2r0 

kJ  ro3  =  Aj  R  K,/2  ir  b 


If »  kj  =  Kj,  rQj  can  also  be  estimated. 


(55) 
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APPENDIX  B:  ANALYSIS  OF  THE  CONICAL  SPIRAL  GROOVE  GAS  BEARING 


The  spiral  groove  bearing  was  tirst  analyzed  by 
RTF  Whipple,  and  his  basic  derivation  was  later  applied 
by  Fortescue  and  Whitley  and  Williams  to  obtain  circular 
thrust  bearing  design  data  Briefly,  Whipple's  derivation 
assumed  an  incompressible  fluid,  negligible  inertia  effects, 
and  a  film  thickness  that  was  small  compared  to  other 
bearing  dimensions.  He  developed  his  equations  for  a 
linear  model  and  assumed  that  adapting  them  to  a  circular 
thrust  hearing  would  result  in  negligible  error  provided  the 
radii  were  large  compared  to  the  width  of  the  bearing  area. 
The  boundary  conditions,  transformed  to  coordinates  along 
and  perpendicular  to  the  grooves,  were 


1  (.  ontinuity  of  tlux  across  the 

groove-ridge  boundary 


,  3  dI’l 
1 


t>u  L  ( li  |  -  h)  cos  0 


(1) 


Periodicity  ol  the  pressure  distri¬ 
bution  along  the  direction  ot  motion 

DP |  DP, 

a  ■  — —  cos  0  *  —  sin  0  = 

1  drj 

•a (-—cost?  sin  o) 

V  3$  dr]  / 


C> 


Figure  10  Inclined  (irt  ><  n  <  Hearing  (icoinctr  i 


3.  The  mean  radial  tlux  across  the 
seal  belt  is  equal  to  the  sum  of  the 
(luxes  in  the  grooves  and  ridges 

i  _p_  ri[  h3dp 
(a,+a)l  2nJ  d> 


(i^T.  it  h  "i  (Iime 

-|jcoSs)*«hJ(|i'sm#|Ea»4#)J 


(3) 


When  these  equations  are  solved  simultaneously  and  trans¬ 
formed  back  to  coordinates  along  and  perpendicular  to  the 
line  of  motion,  then 


1  a.a  sin-0  (Ik -ItV") 

I  +  all3-  J -  > -  l 

alf,  +  aili3  J 


alh 


■r-i1 

(h|-h)  a  |  a  (lf*-h3)  sin  0  cos  0 

U - ! - 

3 


(4) 


all,  +  a  |  It* 


If  two  bearings  are  assembled  back-to-back  so  that  the 
pressure  of  one  opposes  that  of  the  other,  the  net  tlux  is 
zero,  and 


dy 


=  0 


Furthermore,  if  a,  =  a,  this  can  be  written 


LI  (hj-h)(h|-li3)  sin  0  cos0 
(If*  +  h3)-  -sin-0  (li*  •  h3)- 


(5) 


where  h ,  =  li  t  A 

And  since  h.  It,  and  0  arc  independent  of  y,  the  pressure 
above  ambient  at  an  elemental  conical  strip  is 
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t>*i  U  (hj-hl(lylr)  sm0  cos0  y 


(hj  +  rr)“  -sin*  d/  (hj  -  hj)* 

The  differential  normal  force  acting  on  an  elemental  conical 
surface  (Figure  A-2)  Is 

dF  *  P0  rd  0  dy  (7) 

Developing  this  pressure  distribution  about  a  conical 
surface. 


,  ~  dr 

dy  *  - — 
sin  a 


U  *  rto  (1 

Therefore, 

6pwr  (r  -  r)  (h .  -  h)  (h , 3  +  h3)  sin  0  cos  0 

po(M-r)  = - r ~ "j  3,  -  ' - 3  7 rr 

sin  a  (h  j  +  hJ)“  -  sin*  0  (h  ,3  -  h3)*  1 


The  axial  force  necessary  to  deflect  a  pair  of  conical 
bearings  an  axial  distance 


oa  *  6  sin  a 


r  =  r,  0  x  2  rr 


*.•  f  [ 

«n  a  J  J 

r  -  ru  0  =  0 

•r*(r0-r)  £po  (h  -  6, 0)  -  PQ  (h  +  fi.  0)j  d0dr 


Integrating,  and  assuming  that  6  <  <  h 


(ro"  rl)2(ro:  +  2r0ri  +  3ri:)J(7.®) 

a  h 


where  J(S,0): 


10  cos0y2(2  ♦  y)  [(I  +  y  +  y2)(8*  16  7+  13  y2 


[(2  +  3y  +  3y2  +  y3)2 

+  5  y3  +  y1*)  -y*  sin2  0  (3  ♦  3y  ♦  y2)2] 
-  y2  sin*  0  (3  ♦  3  y  ♦  y2)2]  2 


Figure  II 


A  plot  ot  this  (unction  is  lound  in  Figure  A-3  The  radial 
force  necessary  to  deflect  a  "  -sical  bearings  a  radial 

distance  of  *  6  cos  a 


-2r‘(ro-r)  P0 (  h  f  1  -t  cos  0  ]  .0)  do  dr 


Integrating,  and  assuming  that  e  <  <  1 


3frpu  cot  “a  ,  ,  ■>  . 

- 3 - (ro‘r|t*('o‘tVi  +  3f,‘)J(^) 

n 

(16) 


Figure  12.  Load  Function 
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